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o ■ 

r--s ■ The spin-dependent transport of the electrons tunnehng through a resonant tunnehng structure 

^vj ' with ferromagnetic multi-terminal under dc and ac fields is explored by means of the nonequi- 

, librium Green function technique. A general formulation for the time-dependent current and the 

C . time-averaged current is established. As its application the systems with two and three terminals 

in noncoUinear configurations of the magnetizations under dc and ac biases are investigated, re- 
spectively. The asymmetric factor of the relaxation times for the electrons with different spin in 
0^ ■ the central region is uncovered to bring about various behaviours of the TMR. The present three- 

terminal device is different from that discussed in literature, which is coined as a spin transistor with 
source. The current-amplification effect is found. In addition, the time- dependent spin transport 
'TI2 , for the two-terminal device is studied. It is found that the photonic sidebands provide new channels 

fH ' for the electrons tunneling through the barriers, and give rise to new resonances of the TMR, which 

I ' is called as the photon-asisted spin-dependent tunneling. The asymmetric factor of the relaxation 

5^ , times is observed to lead to additional resonant peaks besides the photon-asisted resonances. 
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I. INTRODUCTION 



Giant magnetoresistance (GMR) effect [l( discovered in Fe/Cr multilayers has motivated much research on the spin- 
dependent transport in hybrid magnetic multilayers in recent years^. Parkin et al.'S] has observed that the interlayer 
O \ exchange coupling oscillates damply with the Cr or Ru spacer layer thickness in Co/Ru, Co/Cr and Fe/Cr superlattice 
systems, which suggests that the quantum-well states be formed in the spacer layer. This is because the electrons 
may meet multi-reflection at the interfaces between the ferromagnets and the spacer layer, leading to the oscillation 
of density of states (DOS) of electrons at the Fermi level, and thereby giving rise to the oscillation of the interlayer 
exchange coupling with the spacer thickness|3,la,0- On the other hand, the magnetic tunnel junction (MTJ) such as 
CN \ the ferromagnet/insulator/ferromagnet (FM/I/FM) structure was found to display the spin- valve effect Jj. The MTJ 
^N ■ with double tunnel barriers such as FM/I/FM/I/FM structure was shown to reveal more interesting behaviors on 
^S^ ' the variation of the tunnel magnetoresistance (TMR) with bias^. Recently, advances in nanolithography and thin- 
^1 film processing make it possible to fabricate very small double tunnel junctions called the single electron transistor 

(<— ^ ' (SET). With a very narrow spacer, the Coulomb interaction is important and may lead to new phenomena such as 
■^ Coulomb blockade. Further investigation about the spin-dependent transport in a ferromagnetic SETlyj was performed 

d intensively in the sequential tunneling regime lOj llj and in the cotunneling regime 12, 13]. In the former case, the 
oscillation of the TMR with a dc bias is observed; while the TMR in the latter case is enhanced in the strong tunneling 
I ■ regime, and the Coulomb blockade region is squeezed by the spin accumulation caused by the cotunneling. It would be 
■ O '■ interesting to consider another limit where the size of the spacer is not so small that the Coulomb interaction is weak, 
^ ■ and the quantum well states can be formed in the spacer. If an ac field is applied to the leads and the central spacer, 
one would expect that the interesting features of the TMR might appear due to new effective tunneling channels from 
the photonic sidebandsQ] induced by the ac field. 

On the other hand, Johnson 15] demonstrated that the multi-terminal system is as important as the two-terminal 
system. In the so-called spin dipolar switching device (i.e. the FM/nonmagnetic metal(NM)/FM structure), the 
;J] ■ current flows from one FM layer and out of the NM layer, while the other FM layer referenced to the NM layer serves 
Cu ' as a voltmeter to give an output voltage which depends on the magnetization configuration of the two ferromagnets 15] . 
Fert and Lee|l6j] interpreted Johnson's experiment by means of the Boltzmann equation, and at almost the same time, 
Hershfield et al.[T3] proposed a set of weak coupling equations to describe Johnson's exp eriment. By invoking the 
framework of the scattering matrix theory (SMT) developed by Landauer and Biittiker[l^ for the mesoscopic devices 
with multi-terminal, Brataas et al. 19, 20] discussed a three-terminal device which is to some extent different from 
Johnson's experimental layout, where the current does not flow out of the NM layer but out of another FM layer, and 
the NM layer serves as a central spacer region, while a third FM layer is introduced to measure the chemical potential 
and no net current flows through it. In other words, the current entering the third FM lead is equal to that flowing out 
of it. It is found that the direction of the magnetization of the third FM lead can affect remarkably the current flowing 
through the first and the second FM leads when the magnetizations of these two leads align antiparallel. In this case, 



the intrinsic spin relaxation time in the central spacer was assumed to be sufficiently shorter than the order of the 
time between successive tunneling events such that the spin accumulation efFect,21j could be negligible. Furthermore, 
Jedema et a/. [22, 123 investigated the spin injection and spin accumulation in an all-metal lateral mesoscopic spin 
valve with multi-terminal. All these works were performed under dc biases. As mentioned before, if an ac electric 
field is applied to the magnetic hybrid junction with multi-terminal, one would expect that different spin-dependent 
transport behaviors of the electrons might appear. 

In this paper, the time-dependent spin-polarized transport of a resonant tunneling structure with ferromagnetic 
multi-terminal will be investigated by means of the nonequilibrium Green function technique (NEGFT). It is presumed 
that every terminal and the central spacer are applied by the ac fields. In our model, although the spin accumulation 
could present, as we shall pay our attention to the current and TMR in response to ac fields, the spin accumulation 
will not be considered, which will be given elsewhere. A general current formalism for such a system with N terminals 
which are made of ferromagnetic materials or normal metals in the presence of an ac electrical field will be developed. 
A three-terminal device, to be coined as a spin transistor with source (STS), in which every terminal is applied by a 
source bias, and there is net current flowing through each of them, will be proposed. The amplification effect of the 
electrical curent by changing the direction of the magnetization of the third FM terminal in a STS will be discussed. 
The photon-assisted spin-dependent tunneling in a two-terminal device under an ac bias voltage will be explored. 

The rest of this paper is organized as follows. In Sec. II, a model is proposed and a general formalism for the 
time-dependent current in a resonant tunneling system with ferromagnetic multi-terminal is derived. In Sec. Ill, the 
systems with two and three terminals in dc steady states are discussed. In Sec. IV, the spin-dependent transport 
under an ac bias for a system with two ferromagnetic terminals in which the magnetic moments are noncollinearly 
aligned is investigated. In Sec. V, a brief summary will be presented. 

II. MODEL HAMILTONIAN AND GENERAL FORMULATION 
A. Hamiltonian and the Uncoupled Green Function 

Consider a resonant tunneling structure with ferromagnetic multi-terminal (RTSFMT) with the Hamiltonian given 



by 



-W — Pleads 1 ^center i J^couplingj \^) 



where Hii^ads = l^a^a (ct = 1,2,3,...) is the Hamiltonian of the ferromagnetic leads with H^ the Hamiltonian of 
the ath lead, Hcenter = J2cr^d{t)dl.d„ is the Hamiltonian of the central spacer with Edit) = Eo + ^d{t) the single- 
particle energy in the central spacer modulated by a time-varying external field where Eq is the single-particle level 
in the absence of the external field, Ad{t) is the energy from the time-varying field applied to the central region, and 
dl^{da) is the creation (annihilation) operator of an electron with spin a in the central region. Here we presume that 
the single-particle ener gy l evel, say Sd(t), which can be viewed as a single-impurity level or a single energy level of 
quantum well states |24l l25l l2a | . is spin-degenerate. This case corresponds to a nonmagnetic island where the numbers 
of electrons with spin up N^i and down Ndi are equal. In the equilibrium (e.g. in the absence of an external bias), the 
net moment at the nonmagnetic island is zero, while in the nonequilibrium (e.g. in the presense of an external field) 
the numbers of electrons with spin up and down have different tunneling rates entering into the island, and also have 
different tunneling rates escaping from the island. Thus, if the spin relaxation time is not very short (for example, 
longer than the tunneling time), the spin accumulation will occur. In the present model, we shall pay attention to 
the spin-dependent transport of the current and the TMR of the structure, and we will not discuss the effect of 
spin accumulation here for clarity J23|. H coupling is the coupling interaction between the leads and the central region, 
defined by 

^coupling — / ^ ^Ta^ [/^ ) 

a 

Htc = Y. [Tk:dithlj.' + h.c], (3) 

where 7^ ^ is the creation operator of an electron in the ath lead with momentum ka and spin a, T^^di^) is the 
element of the tunneling matrix which describes the coupling between the leads and the central region, and depends 



on time, spin and momentum. We assume that all leads are ferromagnetic, and one of them is chosen as the lead 1. 
The magnetization of the lead 1 is presumed to align along the z axis, and the magnetizations of the other leads are 
supposed to deviate the z axis by angles ^2, ^3, • • ■ , etc. The Hamiltonian of the lead 1 can be written as 

^i=E4.WaL«fc- (4) 

ha 

where a^^^^iaka) is the creation (annihilation) operator of an electron of the lead 1 with momentum fc and spin 
tr, e\„(f) — £k(t) — (iMi is the single-particle energy modulated by a time-dependent external field with ek{t) = 
£fc(0) -I- Ai(f), (T = ±1, Ml = ■^gfiBhi with g the Lande factor, /is the Bohr magneton, and hi the molecular field of 
the lead 1. For the noncoUinear configuration of the magnetic moments of the leads, the Hamiltonians of the other 
leads read 

^b = J2^[£k,{t) - (jMb cos Obibl^Jjk.a - Mb sin Obbl^Jjk.a}, (5) 

ki^a 

where 6 = 2, 3, • • • denotes the lead 2, 3, • • • , etc., fct, is the electron momentum in the lead b, bj, ^(bk^a) is the creation 
(annihilation) operator of an electron in the lead b with spin a {a = — cr), £kt{t) = £fc(,(0) -I- Af,(t) with Af,(t) the 
energy from the time-dependent external field applied to the lead b, and Mi, = 29I^Bhb with hb the molecular-field of 
the lead b. In terms of Eq.Q, we have the coupling Hamiltonian between the lead 6 (= 2, 3...) and the central region 

HTb = J2 TO' (t)^Lda' + h.c] . (6) 

kb<Ta' 



By performing the u — v transformation 



bkta = cos —bk,a - cr sin —bk^a, 

It ^b,'f . ^b,t 

bi „ = cos — bl — crsm — ol, -, 



to Hb, we find that Hb becomes 



Hb = Y.sl,{t)bl^bk,a, (7) 

with e^ ^{t) — ek^{0) ~ aMb + A(,(i). Consequently, the coupling Hamiltonian becomes 

HTb = E TO'(0(cos ^-^blj.' - asin ^-^bl-d,,) + h.c.]. (8) 



ktacr' 



To this end, the Hamiltonian of the system is well defined. In a time-dependent case, the current conservation can 
be expressed as 



E-?«w 



e ; = 0, a = 1, 2, 3. 

dt 



where j%{t) is the tunneling current flowing from the ath lead to the central spacer, while i''(i) = edNd{t)/dt is 
the displacement current, with Nd{t) — '^^ (d\drj\ being the occupation number of electrons with the energy level 
ed{t) in the central spacer. Recently, there are a few works towards partitioning the displacement current into 
each terminal J23, l29j| . Likewise, the displacement current in the present case is supposed to be partitioned into 
every terminal in such a way, i.e. J'^{t) = J^a'^ai^): that we have J^ai^ai'^) + ^ai^)] = J2aJa{t) — 0, where 
Ja{t) — Ja{t) + Ja{t) is the total current flowing from the ath lead into the central scattering region, and J^{t) — 
— ^ {[H, Na]) = — ^ {[Hto, No]) with Na the particle occupation number operator in the ath lead. By means of the 



■y tne pa.] 



nonequilibrium Green function technique |3(l. l31| . after some calculations, we obtain 



where 



Rn 



/ J.TT j^n \ 

\ ^ka,d ^ka,d / 



and 6a is the angle between the magnetization direction of the ath lead and the z axis. By using the equation of 
motion of the Green function, and then noting the Langreth theorem |3l). we get 

J^(t) = ^Ke^ /rfiirr.[GS(t,ti)S<(ti,i) + G<(t,ti)S^(ti,t)], (9) 

where G^(i,t') — I ^^^^ ' / ^i^^ ' /\ ) ^^ ^^^ lesser function of the central spacer in spin space, with G^(i,t') 
being the retarded Green function of the central spacer in spin space: 






where G^ '^{t-,t') = ^HT Ida (t) d*^, (t')>\. In Eq.®, S2(ii,t2) (7 =<,>,'r,a) is corresponding to the lesser, 
greater, retarded and advanced self-energy, respectively, which is defined as 

^l{h,h)=J2(Qtdih))^slJti,hmtdit2)), 7=<,>,^,a, (10) 

kc 

where g^ (ii , t2) is the corresponding Green functions of uncoupled leads, with the retarded (advanced) Green function 
given by 

g'-J^fit, t') = gll^^it - t') exp(T* /* di"A„(t")), 

Jf 

and the lesser function of the uncoupled leads given by 

^'="^'^"1 /(£^^Jexp(T*et^(t-t')) j 

•exp(Ti / dt"Aa{t")), 
Jt' 

where / is the Fermi function. When the central region is uncoupled to the leads, the retarded Green function of the 
central region can be obtained exactly: 

55(t,0 = -*^(i-i')cxp(-z / dt"[eo + A<i(i")]V 



B. Self-Energy S and the Line- Width Function T 



From Eg. 1)10(1 . we shall calculate the self-energy. The off-diagonal elements of the coupling matrix describe the 

spin-flip scattering processes, which was discussed in Ref.|33l- To focus on the low- temperature properties of the 

transport, where the available electrons near the Fermi level are dominant in the transport process, it is reasonable to 

ignore the effect of the spin-flip scattering during the tunneling process, and the coupling parameters are presumed 

7 T t \ 
to be independent of momentum. Therefore, we have T^: d — { n rr .To perform the momentum summation 

\ ^ J^ai J 

of the electrons, we may first sum the momenta of electrons in the uncoupled Green functions of leads. Then, for the 

retarded and the advanced Green functions of the uncoupled leads, we will obtain the real part A and the imaginary 

part r. Because A makes the energy levels shifted and can be absorbed into the energy level of the central spacer, we 

only consider the imaginary part F. With this consideration, one may get 



S;(i,0 - -^y'^e--(*-'')r„(0„,e,t,Oe-'-^*"^'"'^°(*"\ 



where the line-width function Fq, is 



with 



iT^ 



ii^ 



(11) 



(12) 



F^^(0„,£,i,t') = 27r^T:^(i)T„T(i') 



2 '^a r/ 
cos £ 

2 ^ 



st^)+sm'^S{e~et^^) 



FJ^(0„,e,i,t') = 2^^T„*^(i)T„x(i') 



sin^fS{e-et,) + cos'fdie-etO 



F^^(0„,£,i,t') = 27r^sin^cos^T„*^(i)T„x(t')['5(£-£LT)+'5(^-^Li)] 



Ff^(0,,£,i,t') = 27r^sin^cos^T:^(t)T„T(t')['5(£-eLT)+^(^-£Li)] 



Then, the lesser self-energy is 



de 



where /^(e) is the Fermi function of the ath lead. As a result, the current can be expressed as 



(13) 



r^t) = --3m [df [^e-^eit'-t)^-^Ifdt"A.it")J,^^ 
h J J 2tt 

{r„(0„, £, i', t) [G< (i, t') + 2/„(£)GS(i, t')] } . 
It is nothing but the current flowing out of the ath terminal. 



(14) 



C. Equation of Motion of the Green Functions 

Let us first specify the forms of the energy shifts caused by application of ac biases to the leads and the central 
spacer. They bear the usual forms: 

A„(t) = T4cosKi), (15) 

Ad{t)^VQCos{LUot), (16) 

where Wq(o) is the frequency of the corresponding external ac bias to the ath lead (the central spacer), and e ^ h ^ 1 
is assumed. 

By applying the nonequilibrium Green function technique, we get the equation of motion of the retarded Green 
function 

{idt - £d(t))G5(i, t') = 6{t - t')I+ J dtiS'-(i, h)G2ti,t'), (17) 

where S''(i,ii) = ^^ SJ^(i,ii) with S^ being defined above. Performing the gauge transformation'SSJ 

G5(i, t') = G;j(i, t') exp(-i f dt"Vo cos(woi")), (18) 

Jf 

and introducing the double-time Fourier transform 

F(w,w') = / dtdt'F{t,t')exp[i{ujt~uj't')], 



„, ,, f duj duj' ^, ,^ r , I /M 

F{t,t')= / —^-F{Lo,Lo')eM-i{Lot - uj't% 

Ea. (|17|l is transformed into 

/di /' ~ ~ 

^S'-(c.,c.")G5(c.",c.'). (19) 

S''(a;, Lj") can be gained by applying the same gauge transformation and the double-time Fourier transform to S''(i, t'). 
Then substituting it into Ea. H19|l . we can in principle get the Green function G2(a;,a;'). To simplify the treatment, 
we will assume without loss of g enerality that the elements of the tunneling matrix are independent of time and T^ 
is taken at the Fermi energy [33J. By changing the summation into an integral, we obtain 



[a; - £0 + ^r(0)]G^(c^, u') = 2^5{lo - J), (20) 



_/ r^(O)(l + P„cos0„) r^^(O)P„sin0„ 



where r(0) ^Y.a'^ci'^) wit^i 

Ta (0) = (^ I ry"(0)P„ sin 0„ " ' T1 (d)Vi - k cos'e,) j ' (^^^ 

with r^(o) = 2npl |T„t|^ r«(o) = 2^pO |r„xl', rf^(o) = 2^pIt*^^t^^, rj^(o) = 2^plT*^{r^^ and pi = p^^ + p^^, 

where Pa]{pa[) is the DOS of spin up (down) subband of the ath ferromagnetic lead, while Pa — {Pa^~ Pai) / {901^+ Pal) 
is the polarization of the ath lead. So the retarded Green function of the central region becomes 

GS(a;,t^') = 27r[c^ - eo + ^r(0)]-i(5(c^ - J), (22) 



where the superscript —1 means the inverse of the matrix. By tranforming G^{uj,uj') back, one gets 

G5(t,i') = / ^e-''^(*"*')e-*.''*' '^*"^ocos(^ot")[^ _ ^„ _ -r(0)]-\ (23) 

J 2tt 2 

G'^{t,t') can be obtained by the Keldysh equation 

G<{t,t') ^ J dhdt2G2t,h)-E<{h,t2)G''At2,t'), (24) 

where S^(ii,t2) = X]q ^q (^i' ^2)- With the presumption of F being real, we obtain the current of the ath terminal 
under an ac external bias: 



+5Re r„(0) J2 fpi^)M^^ t)T0iO)AliLo, i)}, 

/3 



(25) 



where the function A^ of the ath lead is defined as 

A„(w,t)=/ dt'e*"(*-*')e'^*'''*"^°'=°"("°*"'G;,(t,i')- (26) 

J —00 

The displacement current, J^{t) = edNd{t)/dt = ed[Trfj^mG^{t,t)]/dt, can be obtained by 

a 

with 

J^(t) = -e^e-J ^f^{Lo)Tr, [A„(a;, t)r„(0)At,(^, t)] ■ (28) 

The current conservation now becomes X^qI-^qIO + -^qIO] = X]q '^a(^) = Oi where Ja{t) = Ja{t) + >^a(*) is the total 
current of the ath lead. Eqs. (|25|l and H28f) are the main results of this paper. Based on them, we can investigate the 
spin-dependent transport in a hybrid junction system with multi-terminal under either a dc bias or an ac bias. In the 
following, we shall discuss the case in a steady state under a dc bias, and then, we shall consider the spin-dependent 
transport in a system with two and three ferromagnetic terminals whose magnetizations are noncollinearly aligned 
under dc and ac biases. 

III. STEADY STATE UNDER A DC BIAS 
A. Two Ferromagnetic Terminals with Noncollinear Configuration of Magnetizations 

Let us consider the structure whose schematic layout is shown in Fig. 1. We assume that the molecular field of 
the left FM lead is aligned along the z axis which is perpendicular to the current direction, while the magnetization 
direction of the right FM lead deviates the z axis by an angle d. A dc bias is applied between the left and the right 
ferromagnets, which causes a difference of chemical potential of the left and the right lead by /i^ — /z/j = eV. In a 
steady state, there should be no charge accumulation in the central scattering region. So the displacement current is 
zero, and J£ — —Jfi holds. For the convenience of calculation, we then define the current flowing through the system 
as 

J^I{J1~Jr)- (29) 
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FIG. 1; (Color online) The schematic layout of a resonant tunneling system with two ferromagnetic terminals. The energy 
level of the central spacer can be modulated by tuning the gate voltage. The potential profile is shown in the lower panel. The 
magnetizations between the left and the right ferromagnetic leads deviate by an angle 9. 



According to Ea. (|25|l . we need to calculate the function AQ,(e,i). In the steady state, we have 



A„(£,t) 



'2ti uj' -£-%{)+' 



(30) 



where r = t — t'. In addition, we have 



GS(c^) = [..-eo + ^r(0)]-\ 



(31) 



where T — T l + T^. By noting that the pole of G^(w) is on the lower half-plane, while the pole of the integral 
function of A.a{e, t) is on the upper half-plane (w' = £ -I- iO'^), one can make use of the residual theorem and perform 
a contour integral on the upper half-plane, and obtain the function Aq(£) in a steady state: 



A„(£) = G^(£). 



(32) 



It is obvious that the function A is the Fourier transform of the retarded Green function in a steady state. From 
Ea. (P5)l . we arrive at 



2e f 



(33) 



where Te//(w) = lTr^[TLG''^{uj)TRG^{uj) +TRG';^{Lu)TLG'^{uJ)] is the effective transmission coefficient (ETC), and 
from Ea. (|12|l . Tlj^ have forms of 



r[(o)(i + PL) 



r[(o)(i-PL) y 



R — 



Tf{O){l + PRCOS0) Tf^iO)PRSm9 
rf:^{O)PRsm0 rf(0)(l-PflCos( 



with rj^^(O) = r|^^(0) (Recall that T is supposed to be real). 

Before making the numerical calculation, it is better to specify the parameters. First, we assume that the tunnel 
junction under interest is symmetric, say, the two FM terminals are made of the same ferromagnets. Then, we may 
take Pl = Pr = P, and F^qJO) = F^qJO) = T-^q), rf^(O) = Fj^. Defining a parameter which is called as the 
spin asymmetry factor 77 = Fi/Fj, which is similar but different to that in Ref.,lll]. we may write Fi/F = ?7/(l -I- r]), 
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FIG. 2: (Color online) The angle 6 dependence of the ETC (a) and the TMR{e) (b) and (c) for different r). The energy is 
scaled by F, and the parameters are set as P = 0.4, eo = 5, and Ej — 8. 



Fj/r = 1/(1 + rj), (Fii/F)^ = 77/(1 + 77)^ with F = F| + Fj^ as an energy scale hereafter. At low temperature and a 
small bias, the tunnel conductance C can be obtained from Ea. (|33|l : C ~ '^Teff{Ef) with Ef the Fermi energy. The 
tunnel magnetoresistance TMR{6) can be defined as TMR{9) = 1 — C{d)/Cp, where Cp = C{6 = 0) is the tunnel 
conductance when the magnetizations of the two FM leads are aligned parallel. 

The 9 dependences of the ETC and TMR{9) for different 77 are presented in Figs. 2(a) and (b). It is seen that the 
ETC {TMR{9)) increases (decreases) with increasing 9 when r] < 0.6, while it decreases (increases) with increasing 9 
when rj > 0.7. TMR(9) becomes negative when 77 is small. The tunnel conductance and the TMR are remarkably 
affected by the spin asymmetry factor 77 which is the ratio of F| and F|, while Ff (F|) is the tunneling rate of electrons 
with spin up (down) from the left or right electrode to the energy level in the central island. Different rj means different 
tunneling rate for up- and down-spin electrons. The 9 dependence of TMR exhibits the conventional spin- valve effect, 
i.e., the conductance is larger in the antiparallel configuration than that in the parallel configuration for larger rj. 
A small rj (for example 77 < 0.7) implies that F| < T i (i.e. the tunneling probability of electrons from the spin up 
subband to spin up subband is less than that from down to down subbands), leading to an inverse TMR (i.e. the 
conductance in the antiparallel case is larger than that in the parallel case), as shown in Fig. 2(b). When 77 is around 
0.65, the TMR is negligibly small and displays a nonmonotonic behavior, as shown in Fig. 2(c). It suggests that 
the coupling strength of the spin subbands leads to the nonmonotonic change of the TMR. When the polarization 
P = 1, i.e. if the leads are made of half- metals, the TMR{9) will be always positi ve, and TMR{tt) = 1, which 
recovers the perfect spin-valve effect. Here it should be pointed out that Bratkovskv|25J has considered a resonant 
tunnel diode (RTD) with electrons tunneling through a single-impurity level which is similar to the present structure, 
where, by invoking a different method, the 9- and bias-dependences of the conductance were obtained. However, the 
asymmetry of the tunneling rates Em) was not considered there while it is included in the present case, which gives 
rise to interesting behaviors, as manifested in Fig. 2. 

The effect of the spin polarization P of the leads on the ETC as well as the TMR{9) as a function of 9 is presented 
in Figs. 3(a) and (b). It can be seen that the ETC and the TMR{9) are affected considerably by P. At = 7r/2, there 
is a crossing point on the curve of the ETC, showing that the ETC, thus the tunnel conductance, does not depend on 
the polarization when the magnetic moments of the leads are aligned perpendicular. As 9 < 7r/2, the ETC increases 
with increasing the polarization at a given 9] when 9 > it/2, the ETC decreases with increasing the polarization, as 
shown in Fig. 3(a). This is easy to understand, because Teff{9) ~ 1 -I- P^coa9. When P — 1 and 9 — n, the ETC 
is zero, showing again a spin-valve phenomenon. This is because in this situation the majority (spin-up) subband of 
the left lead is fully filled, while the minority (spin-down) subband is empty, if the spin-flip scattering is ignored in 
the tunneling process, the transport of the minority electrons in the left lead tunneling into the minority subband of 
the right lead is prohibited |34|. At a given 9, TMR{9) increases with increasing the polarization P, as shown in Fig. 
3(b). This behavior can be simply understood by the fact TMR(9) ~ [2P'^/{\ + P^)] sin^ 9/2. 

To be consistent with the common definition of the TMR in literature, we define a quantity TMR = 1 — Cap/Cp, 
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where Cap is the tunnel conductance when the magnetizations of the ferroniagnets are ahgned antiparalleL In fact, 
such a definition is nothing but TMR{9 — tt) — TMR. In Fig. 4, we show the TMR as a function of the incident 
energy of electrons. It can be observed that the TMR reveals different behaviors for different ry. First, for any ry 
there is a peak aX E = Eq, and the curves have two crossing points at E = Ei and E2, suggesting that the TMR 
at E = Eq, El and E2 is independent of 77. This property may be caused by the resonant tunneling. Second, when 
E < El and E > E2, the TMR changes from negative to positive with increasing ?/, and larger rj, larger the TMR; 
when El < E < E2 except E — Eq, larger 77, smaller the TMR. It is obvious that the TMR is symmetrical to the 
axis E — Eq. This is a characteristic of the resonant tunneling structure. 
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FIG. 3: (Color online) The ETC (a) and the TMR{6) (b) as a function of 9 for different polarization P. The parameters are 
taken as ?7 = 1.0, eo = 5, Ef = 8. 



So far, we have discussed the spin-dependent transport of the structure at zero temperature. Let us now look at 
what happens at finite but low temperature. To proceed, we set /i^ — /x/j = eV and /x^ = 0. In the rest of this 
subsection, we follow the treatment of Rudzinski and BarnasJ23: and consider the bias-dependent energy level of the 
central region as Ed = e^ — xeV, where x = dL/{dL + du), < x < 1, dL and du are the thickness of the left and 
the right barrier. This treatment is corresponding to introduce a constant Hartree potential in the central scattering 
region as discussed in Ref . |3ItI| . As the structure under consideration is symmetrical, we have x = 1/2 in the present 
case. The dc bias dependence of the current and the differential conductance is depicted in Figs. 5(a) and (b). (Note 
that throughout this paper any voltage is all scaled by F/e.) It can be observed that the current increases with 
increasing dc bias for different energy level of the central region. When Eq — 0, the current is a monotonic function of 
V; when Eq becomes larger, the curvature of the current alters with V, and the current rises dramatically at V^ = -^ 
to saturation, as displayed in Fig. 5(a). As the energy of the central spacer is presumed to be positive, its role is like 
a barrier. Thus, one may see that with increasing Eq the current is remarkably suppressed. When the external dc bias 
lifts the Fermi level of the left lead to meet with the resonant energy level, the resonance may occur, which causes 
the current rapidly rising around the resonant position. This character can be clearly seen in Fig. 5(b), in which the 
differential conductance exhibits resonant peaks at Eq = \e.V , i.e. Sd = 0, a typical character of the resonant tunneling. 
At low bias, the differential conductance exhibits an ohmic behavior especially for large Eq, i.e., is independent of the 
bias voltage; as V goes higher, the conductance deviates the ohmic behavior, which is more obvious when Eq is small. 
Since the appearance of the resonant tunneling, the overall behavior of the different conductance in this junction 
structure differs dramatically from the classical behavior, as manifested in Fig. 5(b). The TMR as a function of 
the dc bias voltage for different Eq is shown in Fig. 6. It is found that the TMR decreases with increasing dc bias 
at low bias but increases at high bias. This feature is consistent with the experimental observationl36| although the 
system investigated in the experiment is a single-barrier structure. One may note that for eg = 2 or 5 the TMR 
first decreases to a minimum, then slowly rises with the dc bias. Before reaching the minimum, there is a "shoulder" 
structure for the TMR. This behavior is also noted by Wang et al.J55j, in which they attribute it to the quantum 
resonance. In the present case we believe that it may be caused by the coupling between the central region and the 
leads. This coupling makes the sharp energy levels in the central region extended to those with a finite width, leading 
to that the electrons in the central region can relax to the leads. It may be the finite width of the energy levels in 
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the central region that results in the "shoulder" structure. In practice, the energy levels in the central region can be 
tuned by applying a gate voltage Vg. We may assume that the energy levels in the central region can be modulated 
in a form of e^ = £0 — xeV — Vg . The dependence of the current and the TMR on the gate voltage Vg for different 
bias voltages is depicted in Fig. 7, where Eq ~ 3. It can be observed that with increasing gate voltage the current first 
rises rapidly, and goes to the maximum at Vg = Eq, then decreases rapidly towards zero, as shown in Fig. 7(a). For 
a large V, the current shows a hump behavior. It can be understood as follows. When a dc bias V is applied to the 
leads, the Fermi energies of the left and the right leads are split by eV , i.e. Ei — £'? = eV . When the resonant energy 
level £d is in the region E^ < Ed < Ei, the resonant energy level provides a tunneling channel and the resonance 
occurs, leading to that the current increases rapidly; when Sd > Ej, the energy level in the central spacer serves as 
a barrier and the current decreases with increasing Vg. When the dc bias becomes larger, the width of the resonant 
peaks become wider, giving rise to a hump behavior. The gate voltage dependence of the TMR is shown in Fig. 7(b). 
When the gate voltage is small, for instance V — I, the TMR goes to the minimum at Vg = Eq where the resonant 
peaks of the current just appear. Around the minimum of the TMR, there appear two "shoulders" on both sides of 
Vg = eq for different V . When V becomes larger, the "shoulder" becomes flatter. With increasing the external dc 
bias voltage, the minimum of the TMR becomes a round peak while the "shoulders" become two minima. The origin 
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of the shoulder behavior is the same as that explained above. 



B. Three-Terminal Device Under a dc Bias 

Now we turn to discuss a three-terminal device under a dc bias. The multi-terminal device has been extensively 
investigated within a framework of the mesoscopic device. For instance, a three-terminal device was discussed in terms 
of the Landauer-Biittiker scattering matrix theory (SMT) 18], where a third terminal is introduced to measure the 
chemical potential /ia and the current flowing out of it is set to zero, say, j^ = 0. In this way, the chemical potential 
of the third terminal can be expressed in terms of the chemical potentials of the first and the second terminals. When 
electrons enter into the third terminal and then flow out of it into the scattering region, their phases are randomized. 
For a spin-dependent three-terminal device, Johnson[l2 observed the spin bottleneck effect in a spin transistor. In 
a structure such as a FM/NM/FM multilayer, the polarized current is injected from one of the FM layers into the 
NM layer and flows out of it; the other FM layer and the NM layer serve as a voltmeter to give a voltage output 
which depends on the magnetization configurations of the two ferromagnets. Johnson and SilsbeeJ23| proposed the 
nonequilibrium thermodynamic theory to investigate the charge transport and the nonequilibrium magnetization. 
Fert and Leejlg investigated Johnson's spin transistor device and discussed the effect of junction resistance on the 
spin accumulation based on the Boltzmann equation. Brataas et al.jlD, 2(j| extended the SMT to a ferromagnetic 
three-terminal system FM/NM/FM, where the current flows from the first FM terminal into the NM spacer region, 
then into a second FM terminal, but the net current flowing out of the third FM terminal is set to zero. When the 
magnetizations in the first and the second FM terminals are antiparallel, the current flowing through the first and 
second terminals can be varied with the magnetization direction of the third FM terminal through which there is no 
net current flowing, though. When the magnetization in the third terminal M3 deviates by 7r/2 to the first one Mi, 
the current goes to its maximum. Here we propose a three-terminal spintronic device which is somewhat different 
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FIG. 8: (Color online) The schematic layout of the resonant tunneling structure with three ferromagnetic terminals. 



from the one discussed in Refs.[19|,|20j. We suppose that the third FM terminal is applied by an external dc bias too, 
enabling the net current to flow through it, just like a usual semiconductor transistor. We call it a spin transistor 
with source, whose schematic layout is shown in Fig. 8. In the following, we shall investigate the spin-dependent 
electrical transport in such a STS. 

From Eq. 1)25(1 , the current of the ath lead is 






■E 



duj 

2^ 



iUuj) - /^(c.)]rr^[r„G^(a;)r^G2(. 



(34) 



where /3 runs over the other leads besides a. We re-arrange the current flowing in the leads as a vector form: 

~J = Bit, (35) 
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where 



D = fC - Cf =- / — (/2 - h)C2l (/2 - /3)C23 

^ -^ ^^ ^ (/3 - h)Czi ih - f2)C32 



/l 
f = I /2 

,00/3 

with /i(2,3) = /i(2,3)(w), and 

Ci2 Ci3 

C = I C21 C23 

C31 C32 

where the matrix elements are given by Cafj — Tr^[T aG^{uj)T i3G'^{uj)] with a, (3 the indices of the leads. In general, 
the formulas can be extended to the system with N terminals, and the current will be a vector with N dimensions 
where C is a A^ x A^ matrix. When one or some of the N terminals are made of NM, then P — Q and the tunneling 
matrix reduces to a number. In the present assumption, C is a symmetrical matrix, i.e. C12 — C21, C13 = C31, and 

(^23 = C'32. 

In the subsequent discussion, we call the parallel (antiparallel) configuration when the magnetizations of the first 
and second terminals are parallel (antiparallel). We set the angle between the magnetizations of the third terminal 
and the first one is 9, as indicated in Fig. 8. Then, the current flowing through the second terminal is 



tc.p _ e 
^2 - ?i 



/ ^{[/2(^) - h{^)]Tr„[T'^G'^/{u^)T,G^/{u^)] (36) 

+[/2(c.) - /3(^)]Tr.[rrG2^(^)r3G2'^(^)]}, 



where the superscript P means the parallel configuration. In the antiparallel configuration, the current flowing through 
the second terminal is 

J^^^'^ - j|^{[/2(c.)-/i(c.)]rr.[rrG--^^(c.)riG-^^(c.)] (37) 

+ [h{Lo) /3(^)]Tr.[rrG2^^(^)r3G«'^^(^)]}, 

where the superscript AP means the antiparallel configuration. At low temperature and small bias voltage, the 
difference of the conductance in these two configurations can be obtained by 

AC = C^{e, =9)- C^P{9:i - 9) (38) 



^rr.[(G"'^r^G2^(c.) - G^,^^{u^)TrG^,^^{u)){T^+T,)l 



Define the TMR{9:i) as usual 



TMR{9^) = ^p,, _,, (39) 



AC 
C({9^ 0) 
Tr,[{G^^P{Ef)TPG^f{Ef)-G':i''''{Ef)T^PG^2^''{Ef)){T,+T,)] 



Tr^[TPGY{Ef,9^ = Q)T^G--P{Ef,9^ = 0)] 
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We shall consider the STS device at finite temperature. The terminal current is determined by Ja = — e ( N^ 

Negative Ja indicates that the positive charges tunnel through the barrier and enter into the central region (i.e. the 
electrons tunnel through the barrier from the central region to the leads), and the minus means the reduction of the 
positive charges (i.e. the increase of the negative charges) in the ath terminal. When Ja is positive, the things become 
opposite. For simplicity, we set ^i = fi^ and fJ.i — fJ-2 = eV, while choose fi2 = 0. In this case, the electrons flow from 
the first and third terminals to the second terminal, suggesting that the positive charges tunnel through the barrier 
into the central region, and J2 is negative. As the flowing direction of the positive charges is the current direction. 



and noting J{^ + J| + J| = 
the second terminal as 



i.e. Jf + J!^ ~ — J|, we define the current flowing from the first and third terminal into 



e) = -{J! + J', 



duj 



[/i(^) 



J2) 
h{Lo)]{Tr,[{T, 



2h J 2tt' 
+Tr.[r2G5(c.)(ri+r3)G3(a;)]}. 



Based on this equation, we could obtain the current J2,p{AP)(fi^ 
numerically. 



r3)G^,(c.)r2G3(a;)] 



(40) 



at the parallel (antiparallel) configuration 
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FIG. 9: (Color online) The Q dependence of the current (a) and the current ratio K2,p(ap) (b) for different polarizations and 
magnetization configurations (parallel and antiparallel). The parameters are taken as V = 5, eo = 8, ksT — O.IF, rj = 1.0. 



The dependence of ^2,p(ap)(^3 = ^) under different polarizations is shown in Fig. 9(a), where the parameters 
are assumed the same as those in the proceding subsection, the temperature is taken as ksT = O.IF, and the three 
terminals are supposed to be made of the same ferromagnets. In the parallel configuration, the current J2,p{03 = 6) 
decreases with increasing 0. This is because of the spin-valve effect, the current flowing from the third terminal into 
the second terminal decreases with increasing 9, leading to the total current into the second terminal decreases with 
increasing 9. While in the antiparallel case, the current flowing from the third terminal into the second terminal 
increases with increasing 9 because M3 rotates from the antiparallel to parallel conflguration with respect to M2, 
resulting in that the total current of the second terminal increases with 9. We note that the current in the antiparallel 
conflguration is equal to that in the parallel conflguration when 9 = tt. Because in this case, i.e. the magnetizations 
of the first and third terminals are antiparallel, the state of M2 || Mi but M3 antiparallel to M2 is symmetrical to 
the state of M2 || M3 but M2 antiparallel to Mi, the current in the two states are equal, and the TMR would be 
zero. Besides, we find that regardless of the magnetization configuration, the polarization makes the absolute value 
of the current increases except 9 = tt. 
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To illustrate the effect of the tunnel niagnetoresistance, for this three-terminal device we define the current ratio as 



<^2.P{AP)\^1. 



J2,P{AP){(^3 



f-J'. 



2,P(AP)(^3 — 0) 



J2,P{AP)i&3 — 0) 



(41) 



This quantity reflects the relative change of the tunnel current with respect to the change of the magnetization 
configuration, whose role is similar to the TMR. The angle dependence of the current ratio k,2,p(ap){G) for the 
second terminal is shown in Fig. 9(b). It can be noted that K,2.p{0) is negative for the parallel configuration and 
the absolute value of K2,p{9) increases with increasing 6. For the antiparallel configuration, K2^ap{Q) is positive and 
increases with increasing 9. K2^ap{'^) is larger than the absolute value of K2,p(7r). This difference becomes larger 
for a larger polarization. For example, P — 0.8, K2,ap[t^) ~ 173%, K2,p(7r) « —35%. It suggests that the change 
of the current for the antiparallel configuration is larger than for the parallel configuration. So, the current flowing 
out of the second terminal can be enlarged or reduced by changing the relative orientation of the magnetic moments 
of the third terminal. The so-obtaind current amplification effect is considerably large for a larger polarization at 
the antiparallel configuration. For example, P = 0.6, the change of the current ratio can be as high as 55%. Such a 
current amplification effect is quite different from that in a usual semiconductor transistor. The present results show 
that in a STS, by merely changing the magnetization direction of the third terminal the current fiowing through the 
first and second terminals can be amplified or reduced. Besides, we can define another current ratio, TMR^iO), for the 
STS to investigate the change of the current flowing out of the second terminal when the magnetization configuration 
varies from parallel to antiparallel: 



TMR^ie) = 



J2.P{0) 



(42) 



The 6 dependence of the TMR^(6) for different polarizations is shown in Fig. 10. One may see that the TMR decreases 
with increasing 6. As mentioned above, the states of the parallel and antiparallel configuration are symmetrical at 
9 ~ TV, leading to the currents equal, and then giving rise to the TAIR^{t:) zero. When the polarization is larger, the 
TMR becomes larger. 

Here we would like to point out that the present three-terminal device is different from that discussed by Brataas 
et al.[l9|, |20j. In their device, the third terminal is introduced as an inelastic scattering resource, and there is no 
net current flowing through this terminal. In the present STS device, there is net current flowing through the third 
terminal, and the chemical potential of the third terminal is not changed (i.e. its chemical potential cannot be 
expressed in terms of that of the first and second terminals). Owing to the spin-dependent scatterings, we have 
found that the current fiowing out of the second terminal can be enlarged by changing 6 when the magnetization 
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configuration of the first and second terminals is antiparallel. The enlargement is more remarkable with increasing 
polarization. This three-terminal device is similar to the usual semiconductor transistor. In the latter case, the 
current flowing out of the collector is enlarged by varying the voltage bias between the emitter and the base, while 
in the present case the enlargement of the current can be realized by changing the direction of the magnetization of 
the third terminal. This is of course of the quantum-mechanical origin. Another difference between our device and 
that of Brataas et al. is that the current can be modified in a STS by changing 6 at the parallel configuration, but it 
cannot be tuned in the case of Refs.p^l2(]|. 

However, the result presented in Refs.[l9l V2Sj( can be recovered from our analysis. Suppose that there is no net 
current flowing through the third terminal in the STS, like the treatment of Biittiker. This implies that the chemical 
potential of the third terminal /is can be determined by the chemical potentials of the other two terminals. As a 
result, we get 

J (Q a\ ^ f'^fff \ ft mC'i2(Ci3 -I- C23) + Cl3C'23 /.g^ 

Ue, = 6) = ^j -[AH-/2H] ^-^^^ . (43) 

In addition, it is assumed that —V/2 is applied to the first terminal and V/2 to the second terminal. The current 
ratio K2,p(AP)(()) can be obtained in terms of Eg. M3II. and the result is shown in Fig. 11. One may find that our 
result is fairly consistent with that shown in Refs.jl3,l20|, although quite different methods are used. 

IV. TWO-TERMINAL DEVICE UNDER AN AC BIAS 

A. General Formulation 

Let us turn to investigate the spin-dependent transport in the magnetic tunnel junction with two terminals under 
an ac bias, where the noncoUinear configuration of magnetizations of the FM leads is presumed. Under an ac bias, 
the total current of each terminal is a sum of the "terminal current" and the "displacement current" with the ac bias 
determined by Eas. (|15|l and (|16|l . The total current is 

J-\iJL- Jr) = ^[(J£ - rn) + {Jt - 4)]- (44) 

Again, we need to calculate the fmiction AQ,(tj,t). From Ea. (|23|l . we have 

A„(a;,t) = / di' /— e^(— ")(*'*')e^^''^*"[^°^°'^("°*")-^°^°^("°*")lGS(u;") 
J-00 J 27r 

J l^Y- A / ^^ ^^ ' f45l 



mnm'n' 



with JaC"n') = Jrn{^)Jn{^)JnA^)Jn'{^J, whcrc usc has bccu made of e-""'; = J2 Jn.{x)e'""^ with J^{x) 



is the TO-th order Bessel function. In the same way, we can get 

/rliij' Hlo" 
— — GS(^")r„G3(c.') 
ZTT ZTT 



m— — 00 



(46) 



^ — i[{m — n)uJa-\-{rn' — n )(jJo]i 



(UJ — LU" + 77lUJa — n'uJo + iO+)(w — LU' + ULOa — m'tOQ — «0+) 

The time-average of a time-dependent physical quantity F{t) is defined by 

fT/2 
-T/2 



{Fit)) = hm i / Fit)dt. (47) 



The time-averaged displacement current can be in principle obtained by means of Eqs.J^HJ and l|46(l . Although the 
integral on time can be nonzero only when the term in the square bracket of the exponent in Ea. l46|) is zero, to get 
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the displacement current a differential of Ea. H46(l with respect to t should be made, leading to that the time-averaged 
displacement current is zero. This is also mentioned in Ref.|3l|. If F{t) is a periodic function of time, its integral over 
time should be finite. When T tends to infinity, the time-averaged F{t) will be zero. When (n' — m')ajo — {m — n)u)a 
holds, the first term of the time-averaged ath terminal current in Eq.jSBJ is 



and the second term is 

rainn'n' (3 

rr.[r„GS(c.")r0G2(c.')] 



(w" — cij — muja + n'uQ — iO+)(Li; — u' -i- nuja — rri'ujo — iO+) 

The time-averaged terminal current is zero for other cases. To get {Ja{t)) i we shall perform the integral over us" 
first. By noting that the pole of G^{uj") is on the lower-half plane, while the pole of the kernel function is on the 
upper-half plane lu" = lu + muja — n'ujQ + iO"*", we may perform a contour integral along the upper-half plane, and get 

mnm'n' 

where G5(a;,eo (m, n')) = [lo — £Q(m,n') -f- fr]"-^, and eQ{m,n') — bq — muja + n'tOQ. To evaluate {Ja{t)) , one 
should first perform the integral over uj" , then make a contour integral on the lower- half plane, and get 

{Jait)f^ = -JSRe Y. ^&')/^//3Hrr4r„GS(c^,e(^(m,n')) (51) 

mnm' 7i' [3 

•r/3G3(^,e^(n,m'))}, 

where G^(a;, Eg (n, m')) = [u; — Eg (n, m') — |r]~-^. As (n' — m')ujQ = (m — n)LOa, leading to Eg (m, n') = Eg (n, m'), we 
write £g(m,n') as Eg. After some algebra, the time-averaged current of the ath terminal can be obtained 

(^"(^)) = i E J^TrA[JcCA')U^)'raG2io,e^)rG'^A^,e^)] (52) 

mnm'7i' 



-[2^^Mw)/M^)raG^(^,E^)r/3GS(c^,Eg)]}, 
/3 

For the time-averaged displacement current is zero, by noting (J) = \{{J1{t)) — (^^(i))), one may obtain 
(Ja(i)) - ^ H /'^{Ji(r„')/L(c.')(Tr.[rLG^(c.,Eg^)r^G^(c^,E^)] 

Tnnrn'n' 

+Tr.[r«GS(a;,E^)riGS(^,Eg^)]) 

-MZ'n')fB{^){Tr,[^BG'^d{^.4)^LG':,{u^r4)] 
-|-Tr.[riGS(^,Eg«)r^GS(^,Eg«)])}. 

Define M = TlG-'^TbG^ then M"^ = (G3)^(rH)^(G5)'^(rz,)^. Because T^, r^, G^ and G;^ are all symmetrical 
matrices, we find M-^ — G'l^T nG'^^T l . With this observation, we get 

(^) = J E /£{^L(r„')/L('^)rr.[rLG5(a;,Eo^)rKG3(a;,E^)] (53) 

-^fl(™"n')/fl('^)rr.[rfiGS(c.,Eg«)riGS(u;,Eg«)]}. 

This is the time-averaged current of the system with two FM terminals whose magnetizations are arranged in a 
noncoUinear configuration under an ac bias, which can be used to study the spin-dependent transport of the spin- 
polarized electrons in the presence of an ac field. 
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B. Photon- Assisted Spin-Dependent Resonant Tunneling 
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FIG. 12: (Color online) The ETC (solid lines) and TMR (dash lines) vary with the energy of the central spacer. Where 
a = Vb/a;o, ujo ~ 0.3, P = 0.6 and Ef — 5. (a) a — lA and -q — 10; (b) a — 2 and rj — 10; (c) q = 3 and ri — 10; (d) a = 3 and 
rj = 1. The ETCs are calculated at e = 7r/3. 



In this subsection, we shall consider a particular situation in which only a dc bias is applied between the two 
terminals while an ac bias voltage is applied to the central spacer region such as Ad{t) = Vbcos(a;oi)- In this case. 



JL{m"n') = •^m'(^) ^ JR&n')^ "^' = ""' ^^d Eq = Sq = ffo + mLUQ = e'q. Thc timc-averaged current becomes 



{J) = 






(54) 



where J™(^ 



is the ?77,-th order Bessel function. The effective transmission coefficient is given by T^jf — 
5 X^m '^™(^)^'''T[rLGJJ(ijj,eg)rflG2(cj,eg)]. At low temperature and under a small dc bias, thc integral can be 
performed by noting limT^o,Vdc^o(/L(^) — /fl(^))/eVdc — S{lu — Ef). As a result, we can get thc tunnel conductance 
G = ^TeffiEf) as well as the TMR = {Gp - Gap)/Gp. 

Based on these equations, the numerical calculation can be carried out. The parameters are taken as those in Sec. 
IIIA. The ETC and the TMR as a function of eq is shown in Figs. 12(a)-(d). It is seen that the ETC and the TMR 
exhibit resonant peaks at certain values of Eq, and the curves are symmetric to the axis Eq = Ef. These resonant peaks 
occur at the positions of the photonic sidebands characterized by ni — 0, ±1, ±2, ... which correspond to the shifts 
of the energy levels caused by the ac field in the central spacer. We have found that the resonant peaks of the ETC 
become sharper and higher as the orientations of the magnetizations of the left and the right ferromagnetic leads vary 
from antiparallel to parallel. With increasing 9, the amplitudes of the resonant peaks are noticed to become small. 
In addition, we have observed that the resonant magnitude of the ETC is sensitive to the ratio of the strength and 
the frequency of the ac field, say, a = Vq/luo, which appears as an argument of the Bessel function. With increasing 
a, the main resonant peak of the ETC at Eq = Ef is suppressed, while more resonant peaks with large |r?T,| appear, 
as manifested in Figs. 12(a), (b) and (d). The resonant peaks of the TMR does not almost change with a. The Eq 
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dependence of the ETC and the TMR for different rj is presented in Figs. 12(c) and (d). One may observe that a small 
r] suppresses the resonant peaks while broadens the peaks of the ETC and the TMR. With increasing ij, the resonant 
peaks become more and obvious. In accordance with the definition, a larger 77 means that it is more difficult for the 
electrons with spin down to tunnel through the barrier and then into the leads, while it is easy for the electrons with 
spin up. So, the resonant peaks of the ETC become even sharper at smaller rj in the antiparallcl configuration, while 
the resonant peaks of the TMR become more obvious. One may note that there are several other peaks of the TMR 
besides the main resonance and the photonic sidebands as shown in Fig. 12(c), which are not caused by the additional 
resonant energy levels of the central region, but by the asymmetric factor of 77 on the spin-dependent transmission. 




Photonic sidebands 
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FIG. 13; (Color online) The schematic illustration of the photon-assisted tunneling. A very small dc bias is applied to the leads 
and an ac field is applied to the central spacer. The quantum well states in the central spacer will be modulated by the ac 
field. The dashed lines indicate the photonic sidebands besides the main resonance which is indicated by the solid line. In this 
illustration, only is the antiparallel configuration of the magnetizations of the left and the right ferromagnetic leads shown. 

The tunneling process assisted by photons is schematically illustrated in Fig. 13. It can be seen that the photonic 
sidebands provide new channels for the electrons in the left ferromagnet to tunnel through the central spacer and into 
the right ferromagnet. Without the spin-flip scattering, the electrons with spin up (down) in the left will be accepted 
by the spin- up (-down) subband in the right. Because of the splitting of the energy bands of electrons with different 
spins, the spin-valve effect will appear. As £0 increases, the photonic sidebands with negative m, the main resonant 
energy level and the photonic sidebands with positive m pass through the Fermi level of the leads one by one. It is 
these additional channels from the photonic sidebands to make the ETC and the TMR be resonant. 
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FIG. 14: (Color online) The a dependence of the ETC (a) and the TMR (b) for different energy levels in the central spacer. 
The parameters are taken as ujo = 0.3, P = 0.6, rj = 10, and Ef = 5, where 6 = n/3 for (a). 



The a dependence of the ETC and the TMR for different Eq is shown in Figs. 14(a) and (b). One may observe that 
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(i) the ETC oscillates damply with increasing a, and a larger a may suppress the transmission; (ii) the TMR oscillates 
with increasing a; (iii) the peaks are sharper and more obvious when Eq closes to the main resonance and the photonic 
sidebands; (iv) when Eq leaves the main resonance energy to a lower value, the first resonant peak of the ETC and 
the TMR moves towards the direction of larger a's. These characters are the consequences of the photonic sidebands 
induced by the ac field. Here we should point out that the scattering mechanisms of electrons tunneling process 
assisted by phonon and by photon are fundamentally different. Interestingly enough, however, the features of photon- 
assisted and optical phonon-assisted tunneling current and conductance are similar, because both tunneling processes 
are assisted by sidebands induced by photons or optical phonons. In the conventional semiconductor tunnel diodes, 
the tunneling transitions are dominant by emission of phonon at low temperature (see Ref. 38J), while the present 
tunneling transitions are contributed by emission or absorption photon quantum. This photon-assisted tunneling in 
superconductor tunnel junctions was discussed in Ref. 39]. The resonant level is modulated not just by a displacement 
by emitting or absorbing photons, but is modulated in terms of a probability amplitude which is characterized by 
the square of the nth Bessel function J„(Vo/wo) for each level to be displaced in energy by ncoo- These probabiHty 
amplitudes contain the information about the quantum interference. The tunneling probability of electrons via these 
virtual energy levels may be modulated by these quantum interferences. In addition, the photon- assisted tunneling 
can be controlled by tuning an external ac field, but the phonon-assisted tunneling is primarily determined by the 
material. In the present case, the photon-assisted spin-dependent resonant tunneling is considered and the photon- 
assisted TMR is observed. This effect may be used to give a tunable TMR, which may be useful for future application. 
Recenth^ we have noted that a similar phonon-assisted tunneling through a molecular quantum dot was considered 
in Ref.jlg. 

V. SUMMARY 

We have discussed the spin-dependent transport in a resonant tunneling structure with ferromagnetic multi-terminal 
under dc and ac fields by means of the nonequilibrium Green function technique. The general formulation of the time- 
dependent spintronic transport in this structure in the presence of ac and dc fields has been established, which might 
offer a fundamental basis for further discussions on the spin-dependent transport in such spintronic devices in an ac 
field. 

First, we have considered the resonant system with two FM terminals in which the magnetizations of the leads are 
noncoUinearly aligned under a dc bias voltage, and have found that for small ry's where rj characterizes the asymmetry 
of the relaxation times of the electrons with different spin in the central region, the TMR is negative, and whose 
absolute magnitude increases with increasing 9 (the relative angle between the magnetizations of both leads) ; when 
T] becomes larger, the TMR will be positive and increases with increasing 6. We have also investigated the TMR as a 
function of energy, the dc bias voltage and the gate voltage for different polarizations and energy levels of the central 
scattering region, respectively. The results are diverse, which manifest the effect of the resonant energy level on the 
TMR and the tunneling current, as discussed in detail in Sec. II. 

Second, we have considered a three-terminal device under a dc bias, which is different from the standard three- 
terminal device discussed elsewhereRefs. p^ l20j . In the present device, each terminal is applied with a source bias, 
suggesting that there is net current flowing through every terminal. Regardless of the spin-flip scattering, the electrons 
in up-spin subbands in one terminal will be accepted by the up-spin subbands in another terminal. When the 
magnetizations of the first and second terminals are antiparallel, the change of 9 which is the relative angle between 
the magnetizations of the first and third terminals may considerably change the current fiowing out of the second 
terminal. It turns out that tuning 9 the current can be remarkably enhanced, giving rise to the so-called magnetization 
configuration- induced enhancement of the current or the current ratio K2.p{AP)i(^)- It has been found that the 
TMR-i{9), which is defined in this case by the current change ratio when the configuration of the magnetizations of 
the first and second terminals is turned from parallel to antiparallel, decreases with increasing 9. 

Third, as an application of our general formulation we have investigated the time-dependent spintronic transport 
in a system with two terminals under an ac field. We have uncovered that, when a very small dc bias is applied to 
the two terminals but an ac field is applied to the central scattering region, the photonic sidebands are formed in the 
central region. When the sidebands meet with the Fermi energy of the terminals, the photon- assisted tunneling will 
occur. It manifests as the resonant peaks in the ETC versus eg- The TMR exhibits many resonant peaks at the main 
resonant level and the photonic sidebands. It has been found that the asymmetric factor 77 may lead to additional 
peaks besides the photon-assisted resonant peaks. We would like to mention that the present study might open a way 
to control the spin-dependent transport in a spintronic device by applying a time-dependent electrical field. 
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